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^ Abstract 

l/^ ' We obtain an asymptotic expansion for the solution of the Cauchy problem for 

the Korteweg-de Vries (KdV) equation 

ut + 6uux + e'^Uxxx ^ 0: M(x,t = 0, e) = Mo(x), 

for e small, near the point of gradient catastrophe {xc,tc) for the solution of the 
dispersionless equation Ut + duu^ = 0. The sub- leading term in this expansion is 
' described by the smooth solution of a fourth order ODE, which is a higher order 

analogue to the Painleve I equation. This is in accordance with a conjecture of 
Dubrovin, suggesting that this is a universal phenomenon for any Hamiltonian 
perturbation of a hyperbolic equation. Using the Deift/Zhou steepest descent 
^ ' method applied on the Riemann-Hilbert problem for the KdV equation, we are 

"sO ! able to prove the asymptotic expansion rigorously in a double scaling limit. 

(N 
cn 

! 1 Introduction 



It is well-known that the solution of the Cauchy problem for the Hopf equation 

Ut + Guu^ = 0, = 0) = no(x), x G M, t G M+, (1.1) 



- T— I ■ 

^ ' reaches a point of gradient catastrophe in a finite time. The solutions of the dissipative 



and conservative regularizations of the above hyperbohc equation display a considerably 
different behavior. Equation (jl.ip admits a Hamiltonian structure 

Ut + {u{x),Ho} = ut + ^^^^ = 0' 
with Hamiltonian and Poisson bracket given by 

Ho= u^dx, {u{x),u{y)} = 5'{x - y), 



respectively. The conservative regularization of (jl.ip can be obtained by perturbing 
the Hamiltonian Hq to the form H = Hq + eHi + e^H2 . . . [23j (see also [TOt 119]). 
All the Hamiltonian perturbations up to the order e'^ of the hyperbolic equation (jl.ip 
have been classified in |17l l30] and the Hamiltonian is equal to 



dx. 
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where c{u), p{u) are two arbitrary functions of one variable. The equation ut + 
{u{x), H^} = takes the form 

Ut + e-UUa,. + — [2cUxxx + 'ic'UxUxx + c'ul] + 6^ [2p U^xxxx 

(1.2) 

+2p' (bux^u^^^ + SuxUxxxx) +p"{7uxul^ + Guluxxx) + 2p"'uluj;x] = 0, 

where the prime denotes the derivative with respect to u. For c(n) = 12, p{u) = 0, one 
obtains the Korteweg - de Vries (KdV) 

Ut + 6uux + e^Uxxx = 0, (1.3) 

and for c{u) = 48n and p{u) = 2u, one has the Camassa-Holm equation [7] up to order 
e^. For generic choices of the functions c{u), p{u) equation (jl.2p is apparently not an 
integrable PDF, but however it admits an infinite family of commuting Hamiltonians 
up to order O(e^). 

The case of small dissipative perturbations of one-component hyperbolic equations 
has been well studied and understood (see [4j and references therein), while the behavior 
of solutions to conservative perturbations (jl.2p to the best of our knowledge has not 
been investigated after the point of gradient catastrophe of the unperturbed equation 
except for the KdV case. The solution of the Cauchy problem for KdV in the limit e — > 
0, has been studied in the works of Gurevich and Pitaevskii, [26], Lax and Levermore 
[29j . Venakides [36], and Deift, Venakides and Zhou [lH|T5]. The asymptotic description 
of [Ml [BJ gives in general a good approximation of the KdV solution, but it is less 
satisfactory near the point of gradient catastrophe for the Hopf equation (jl.ip . Before 
this break-up time, solutions to the KdV equation are, in the small dispersion limit 
e — > 0, well approximated by solutions to the Hopf equation. After the time of gradient 
catastrophe, solutions to the Hopf equation cease to be well-defined for all x, while for 
any e > 0, solutions to the KdV equation remain well-defined for all x and t and are 
characterized by an oscillatory region after the gradient catastrophe. In this oscillatory 
region, small dispersion asymptotics for KdV solutions turn out to be elliptic [36]. The 
transition from the asymptotic regime described by (jl.ip to the elliptic regime has not 
been rigorously described yet in the literature. For numerical comparisons we refer to 

m 



This problem has been addressed by Dubrovin in [T7], where he formulated the 
universality conjecture about the behavior of a generic solution to any Hamiltonian 
perturbation of a hyperbolic equation (jl.ip near the point {xc,tc,Uc = u{xc,tc,0)) of 
gradient catastrophe for (jl.ip . Dubrovin argued that, up to shifts, Galilean transfor- 
mations and rescalings, the behavior of the solution of (jl.2p as e ^ near the point 
of gradient catastrophe for the Hopf equation (jl.ip essentially depends neither on the 
choice of the solution nor on the choice of the equation. Moreover, the solution near 
the point (xc, tc, Uc) is given by 

where cq is a constant that depends on the equation, while k depends on the initial 
data, and U = U {X, T) is the unique real smooth solution to the fourth order ODF 



X = TU 



^f^' + ^(U'x + 2UUxx) + ^Uxxxx 



(1.5) 
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which is the second member of the Painleve I hierarchy. In what follows, we will call 
this equation Pf. The relevant solution is characterized by the asymptotic behavior 

U{X,T)=^{<o\X\)^l'^^\<o'^''^T\X\~'^!'^ + 0{\X\-^), asX^icx), (1.6) 



for each fixed T G M. The uniqueness of a smooth solution to (jl.5|) for all X, T E M 
satisfying ()1.6p follows from earlier results [32l [33] , while the existence has been proven 
in[8]. 

The aim of this paper is to prove rigorously that the expansion (jl.4p holds indeed in 
the particular case of the small dispersion limit of the KdV equation near the point of 
gradient catastrophe for the Hopf equation (jl.ip . More precisely the solution u{x,t,e) 
of the KdV equation in the neighborhood of {xc, tc, Uc) has an asymptotic expansion as 
follows, 

I . ^ , ('^(^'^V'\r( x-Xc-Qu^{t-tc) Q{t-tc)\ ,^/4/7\ n v^ 
m(x, t,e) ^Uc + \-pr] U \ J , Y + O e ' • (1-7) 

This expansion holds in the double scaling limit where we let e — > and at the same 
time X Xc and i ^ tc in such a way that 

X — Xc — Guc(t — tc) Q(t — tc) 

hm r-TTT = A , lim - — ^ = T, 

(8A;e6)i/7 ' (4Pe4)i/7 

with X, r € M. The constant k is given by 

k = -f'l'iuc), 

where /_ is the inverse function of the decreasing part of the initial data uq{x), which 
is assumed to be real analytic and with a single negative bump. 

The universality conjecture of Dubrovin should be seen in comparison to the known 
universality results in random matrix theory. For large unitary random matrix ensem- 
bles, local eigenvalue statistics turn out to be, to some extent, independent of the choice 
of the ensemble and independent of the reference point chosen [11^ [T2l I13j . Critical 
break-up times occur when the eigenvalues move from a one-cut regime to a multi-cut 
regime. These transitions can take place in the presence of singular points, of which 
three different types are distinguished [12j . Singular interior points show remarkable 
similarities with the leading edge of the oscillatory region for the KdV equation, while 
singular exterior points should be compared to the trailing edge of the oscillatory re- 
gion. Our focus is on the point of gradient catastrophe, i.e. the break-up point where 
the oscillations start to set in. This situation is comparable to a singular edge point in 
unitary random matrix ensembles. It was conjectured by Bowick and Brezin, and by 
Brezin, Marinari, and Parisi O [6] that local eigenvalue statistics in this regime should 
be given in terms of the Painleve I hierarchy. In [9], it was proven that indeed double 
scaling limits of the local eigenvalue correlation kernel are given in terms of the Lax 
pair for the Pf equation (jl.5p . In addition, an expansion similar to (jl.4p was obtained 
for the recurrence coefficients of orthogonal polynomials related to the relevant random 
matrix ensembles. 

In the setting of PDEs a universality result similar to the KdV case has been con- 
jectured for the semiclassical limit of the focusing nonlinear Schrodinger equation |18j. 
In this case the limiting equations are elliptic and the focusing nonlinear Schrodinger 
equation in the semiclassical limit is considered as a Hamiltonian perturbation of the 
elliptic system. 
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1.1 Statement of result 

Our goal is to find asymptotics as e — > for the solution u{x, t, e) of the KdV equation 

ut + 6uux + e'^Uxxx = 0, u{x, 0, e) = uq{x), 

when X, t are close to the point and time of gradient catastrophe Xc, tc for the Hopf 
equation 

Ut + 6uUx = 0. 

By the method of characteristics, the solution of the Hopf equation takes the form 

u{x,t) =uo{0, x = 6tuo{0 + ^^ 
so that 

u' (^) 
1 + 6tuo(0 

from which one observes that a gradient catastrophe is reached for 
1 



max5giR[-6ii'o(^)] ' 
We impose the following conditions on the initial data uq. 
Assumptions 1.1 

(a) uo{x) decays at ±00 such that 

"+00 

uo{x){l + x'^)dx < 00, (1-8) 

(h) uq{x) is real analytic and has an analytic continuation to the complex plane in 
the domain 

S = {z £ C : \lmz\ < tan6'o|Rez|} U {z £ C : \lmz\ < a} 
where < 6q < Tr/2 and a > 0, 
(c) uq{x) < and has a single local minimum at a certain point xm, with 

u'q{xm) = 0, u'q{xm) > 0, 
Uq is normalized such that uq{xm) = — 1- 

Condition (a) is necessary to apply the inverse scattering transform [3J, while condi- 
tion (b) is requested in order to get some analyticity properties of the reflection and 
transmission coefficients for the scattering problem (see below in Section [2]). Condition 
(c) is imposed in order to have the simplest situation in the study of the semiclassical 
limit of the reflection coefficient of the associated Schrodinger equation. 

Let XM{t) be the x-coordinate where u{x,t) reaches the minimum value —1. Then 
for t < tc we have that 

^-^^'^^ = 6t + /:Kx,t)) <°' forx<XMW, (1.9) 

Ux{x,t) = „^ , , — TTT > 0, forx>XM(i), (1-10) 

6i + f.[u{x,t)) 
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where f± are the inverses of the increasing and decreasing part of the initial data uq(x), 



f±{uo{x))=x, /_(-!, 0) = (-00, xm), /+(-!, 0) = (xm, +00). 

Since f'_ is negative, it follows that there exists a time t = tc for which (jl.9p goes to 
infinity. This happens at the time 

te = -^ max (1.11) 

The above relation shows that the point of gradient catastrophe is characterized also 
by 

f'l{nc)=0. 

Since we have assumed that Uq{xm) 7^ Oi it follows that Xc < XMitc)- Summarizing, 
the point of gradient catastrophe is characterized by (jl.lip and the three equations 

Xc = &t,U, + f_{Uc), 6te + /iK) =0, fL{Uc)=Q. (1.12) 

The point of gradient catastrophe is generic if 

f"{uc) + 0. (1.13) 

Our main result is the following. 

Theorem 1.2 Let uq{x) he initial data for the Cauchy problem of the KdV equation 
satisfying the conditions described in Assumptions and satisfying the genericity 
assumption U.13\) . Write Uc = u{xc, tc,0), with Xc and tc the point and time of gradient 
catastrophe given by \1.12\) and Now we take a double scaling limit where we 

let e ^ and at the same time we let x — > Xc and t ^ tc in such a way that, for some 
X,T eR, 

li- ^ " " ''^ = X, lim ^^f^ = T, (1.14) 

where 

k = -f'!!{uc), 

and f- is the inverse function of the decreasing part of the initial data uo{x). In this 
double scaling limit the solution u{x,t,e) of the KdV equation ^.3\) has the following 
expansion, 

t, .)=u,+('-f] u ( .°:^!'-'-' . + o f.^/^ . (1.15) 



{8ke^)7 (4fc3e4 



Here U {X, T) is the unique real pole-free solution to the Pf equation il.5]) satisfying 
the asymptotic condition l[1.6\) . 

Remark 1.3 The expansion (ll.lSp is very similar to the one obtained in [9] for the 
recurrence coefficients of orthogonal polynomials with respect to a weight e~"'^^^^ on 
the real line, in a double scaling limit where the potential V tends, in a double scaling 
limit, to a critical potential with a singular edge point. 

It should also be noted that the error term in (ll.lSp is of order e^^''. It will follow 
in a nontrivial way from our analysis that the term of order e^/'', which a priori seems 
to be present when proving Theorem II. 2^ vanishes. 
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Remark 1.4 If we would consider initial data that do not satisfy the generic condition 
(jl.l3p . our result is not valid any longer. It is likely that the role of the solution U 
would then be taken over by a smooth solution to a higher member of the Painleve I 
hierarchy. 

The proof of our result goes via the Riemann-Hilbert (RH) approach. The starting 
point of our analysis will be the RH problem for the KdV equation, developed in [3l [35] 
using inverse scattering. The Deift/Zhou steepest descent method [16] has shown to 
be a powerful tool in order to obtain asymptotics for solutions of RH problems. The 
strategy of this method is to simplify the RH problem in several steps by applying 
some invertible transformations to it. At the end this leads to a RH problem for 
which asymptotics can be easily found. The Deift/Zhou steepest descent method has 
not only been fruitful in the field of integrable systems, it lead also to universality 
results in random matrix theory [IH [T2l [T3] and to remarkable combinatorial results 
such as in [2j. For the KdV equation, a steepest descent analysis was carried out by 
Deift, Venakides, and Zhou |14i I15j . We will follow the main lines of their approach, 
but however with some important modifications which are necessary to perform the 
analysis near the point of gradient catastrophe. The first modification is to perform 
the steepest descent analysis of the RH problem associated to the initial data uq (x) and 
to the corresponding reflection coefficient itself, while in [14^115] the reflection coefficient 
was identified with its WKB approximation before carrying out the RH analysis. Also in 
the works by Lax and Levermore |29j the analysis was performed on some approximate 
initial data UQ{x,e) for which tto(x,e) uq{x) as e ^ 0, such that the reflection 
coefficient identifies with its WKB approximation. The second modification concerns 
the so-called ^-function, which we need to modify in order to have a RH problem that 
behaves smoothly in the double scaling limit. The third and probably most essential 
new feature is the construction of a local parametrix built out of the \I'-functions for 
the equation. 

1.2 Outline for the rest of the paper 

In Section [2l we give a short overview of the inverse scattering approach in order 
to arrive at the RH problem for the KdV equation. We also recall some previously 
known asymptotic results on the reflection and transmission coefficients, which will be 
necessary in Section [3] when performing the Deift/Zhou steepest descent analysis of 
the RH problem. In the asymptotic analysis of the RH problem, we will construct a 
^-function which is slightly modified compared to the one used in [14^ I15j . After the 
opening of the lens, the crucial part of Section [3] consists of the construction of a local 
parametrix near the critical point Uc- Here we will use a model RH problem associated 
to the Pj equation. Accurate matching of the local parametrix with the parametrix 
in the outside region will provide asymptotics for the solution of the RH problem. In 
Section [4] finally, we collect the asymptotic results obtained in the previous section to 
prove Theorem [121 

2 Inverse scattering transform 
2.1 Construction of the RH problem 

In the first part of this section, we recall briefiy the construction of the RH problem 
associated to the KdV equation, as it was done in e.g. [H [35]. The RH problem 
is constructed by inverse scattering and will be the starting point of our asymptotic 
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analysis in the next section. Recall that we consider initial data satisfying Assumptions 
0(a)-(c). 

The initial value problem for the KdV equation 

ut + 6uu::, + e^uxxx = 0, uix,0,€) = uq{x), (2.1) 

can be solved by inverse scattering transform [22] . Introducing the operators L and A, 
depending on x and also on t through u = u{x,t), 

n (f . . ') ( d d 

L = e^—, + u, A = i€^—r + 3(u— + —u 
dx'^ dx'^ \ dx dx 

the KdV equation can be written in the Lax form [28j 

L = [L,Al (2.2) 

where L = — is the operator of multiplication by ut- The Schrodinger equation with 
dt 

potential u 

can be re-written as the first order differential equation 
DF = JF + QF, D = e^, 

There are two fundamental 2x2 matrix solutions ^ = ^'(A; x, t, e) and $ = ^{X; x, t, e) 
to the Schrodinger equation ()2.3p . with the following normalizations as x — > ±00, 



*(A; X, t, €)e-^-^> ^'"3 _ A(A), as a; ^ +00, 
$(A;x,t,e)e^(-^)''"^'"3 ~ A(A), as x ^ -00, 



where 



A(A) = (_^(_\),/2 .(-A)V2)> -3=(; _"^). (2.5) 

Here we take the principal branch of the square root, such that (— A)^/^ is analytic in 
C \ [0, +00) and positive for A < 0. We denote by ^'i, ^'2, 'I'l, and $2 the columns of 
the matrices ^ and $ respectively. It is known that the vectors ^1 and $2 are analytic 
for ImA > and continuous for ImA > 0, and that ^2 and $1 are analytic for ImA < 
and continuous for ImA < [31j . 

For analytic initial data satisfying Assumption ll.ll fb). ^'(A;x,t,e) and ^{X;x,t,e) 
are also analytic for (x, A) £ S x Hq^ [20], where S is defined in Assumption 11.1( b). and 

Ue^ ={AgC\{0}, -200 + vr < argA < 200 + vr}. 



Denoting ^'(x. A)* = ^(x. A), it follows from the reality of ^0(2;) that 
$2 = ^'i, ^>2 = ^>t. 
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Since \I/ and ^ form two fundamental bases of solutions to the Schrodinger equation, 
they are related to each other by a constant matrix (independent of x) for A G Ilgp, 



^{X;x,t,e) = $(A;x,t,e) 



a(A;t, e) b{X;t,e) 
b{X; t, e) a(A; t,e) J 



(2.6) 



The components o(A;t, e) and b{X;t,e) (and also a*, b*) are holomorphic with respect 
to A S Hog- Since det ^' = det $, we have that 



aa* - bb* = 1, 
and, for real A < 
\af-\bf = l, 

which shows that a ^ for real A < 0. Therefore we can divide by a obtaining 



(2.7) 
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~ < 



1 l(_A)i/2:, 



-e e 



as X ^ — oo, 
as X — > +00. 



(2i 



a 



The quantities 



r(A;i,e) := 



KA;t,e) 
a(A;t,e)' 



t(A;t,e) := 



a(A;t,e)' 



are called reflection and transmission coefficients (from the left) for the potential 
u{x,t,e). For an arbitrary dependence of u{x,t,e) on t, it is not possible in general to 
find the time-dependence of a and b. If however u{x, t, e) evolves according to the KdV 
equation, the Gardner-Greene-Kruskal-Miura equations (22] 



^ = 0, ^ = ^(-A)3/25, 
dt ' dt ' 



(2.9) 



hold so that the reflection coefficient evolves according to 

r(A;M)=r(A;0,e)e?(-^)^^^*. 

In what follows we write r(A; e) = r(A; 0, e) for the reflection coefficient at time t = 0. 
For A —oo we have 131 



lim r(A; e) = 0. 



Remark 2.1 For analytic initial data satisfying Assumptions [TTTl the reflection coef- 
flcient r(A; e) and the transmission coefficient i(A; e) are meromorphic functions in the 
sector Ilgg. The possible poles occur at the zeros of a(A;e). For initial data satisfying 
Assumption ll.lt the zeros of a(A; e) lie in the sector [21] 

vr < arg A < vr + 29o . 

The time evolution of the scattering data now leads us to the inverse scattering 
problem to recover the KdV solution u{x, t, e) from the reflection coefficient at time t. 
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This brings us to the RH boundary value problem for the initial value problem of the 
KdV equation. The 2x2 matrix-valued function M defined by 



M{X;x,t,€) 



V a{X;t,e) J 



a{X;t,e) 



(2.10) 



where and $j denote the j-th column of "3^ and $, satisfies the following RH condi- 
tions [3135]. 

RH problem for M 

(a) M(A; x, t, e) is analytic for A G C\M, 

(b) M has continuous boundary values M_|_(A) and M_(A) when approaching A € 
M \ {0} from above and below, and 

/ -1 (\ \ 2ia{X:x,t)/e\ 

M+(A) = M_(A) ' a(A:.,t)/e 7 I for A < 0, 



f(A;e)e-2i"(^;^^*)/e i_|r(A;e)|2 
M+(A) = M_(A)c7i, ai = for A > 0, 



(c) M(A;x,t,e) ~ I |, for A ^ oo. 

liV— A — iV— A> 



Here a is given by 

K(A;x,t) =4t(-A)3/2 + 2;(-A)i/2_ 



The solution of the KdV equation can be recovered from the RH problem by the 
following formula, see e.g. [3], 

u{x,t,e) = -2iedxMl^{x,t,e), (2.11) 
/Vfl (r f f) 

where Mn (A; x, i, e) = 1 + ^ + 0(1/A) as A ^ oo. 

2.2 Semiclassical limit of the reflection and transmission coefficients 

In order to study the limit as e — > of the Riemann-Hilbert problem, it is necessary to 
obtain semiclassical asymptotics for the reflection coefficient r(A;e). 

Under the validity of Assumptions 11.11 the WKB approximation of the reflection 
coefficient as e — > is the following, see [3^ [20]. 

(i) For any positive constant 5, we have the following WKB asymptotics for — 1-1-5 < 
A < 0, 

r{X;e) = ie-TPW^l + eh^{X;e)), 

t{X; e) = e-"W/^e^(''(^)-'^(^))/^(l + e/i2(A; e)). 
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Here p{X) and p{X) are defined by 

, r^-W , , 1 f (u)du 

p{X) = x_{X)V^+ / yuo{x)-X-V^]dx = - / ^-U^, (2.13) 

1 fU+{u)du 



p(A) = x+(A)^- / [Vuo(x)-A-^]dx = - / ^^±^,(2.14) 

where x_(A) < x+(A) are the solutions of the equation uo{x±{X)) = X and /^i 
are the inverse function of the decreasing and increasing part of the initial data 
respectively. The function r is defined as 



r(A) = J yJX-UQ{x)dx = J VX-uf'^{u)du+ J VX-uf'_{u)du. (2.15) 

The functions /ii and /12 are classical analytic symbols of nonnegative order 
The expansion (j2.12p can be extended to a strip of the complex plane of the form 

n+ = ne„ n {-1 + - < Re (z) < 0} n {0 < Im (z) < a}, (2.16) 

where < a < cr, where a has been defined in Assumptions fTTTl The constant a 
has to be chosen sufficiently small so that for A G 11^ there exist only two simple 
turning points X-|-(A) G 5 solving the equation uq{x±{X)) = X. 

(ii) For any fixed 6 > 0, there is a constant C2 > such that for X < —1 — 6, we have 
that 

r(A;e) = 0(e-^), as e ^ 0. (2.17) 

(iii) Near —1, for sufficiently small 5 > 0, we have the following asymptotics for 
A E D+{-l,6) = {A e C : |A + 1| < 5, ImA > 0}, 

. exp[-2ip(A)/g] , ^. „ 
r(A;e)=z / ir{X)^ i^ + 0{e)), as e ^ 0, (2.18) 



N 
where 
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N{z) = ^L?Z!_e^iog(^/e) (2.19) 
^ ^ r(l/2 + z) ' ^ ^ 

and r(z) is the standard F-function. The functions p(A) and t(A) are the analytic 
continuations to -D_|_(— 1, 6) of the corresponding functions defined on the interval 
-1 < A < 0. 

Remark 2.2 For complex values of A close to —1, the two simple turning points 2;±(A) 
near x = xm are x±{X) XM^s^/X^\^, with = 2/u'q{xm)- In particular for A > — 1, 
these turning points are real, and for A < — 1, they are purely imaginary. Choosing the 
branch cut of p for A < —1, we can, for 6 sufficiently small, extend p to an analytic 
function in D-|_(— 1,(5). 



^ A function h{z\ e) defined in (7 x (0, eo) wfiere U is an open set in C, and eo > 0, is called a classical 
analytic symbol of order m if /i is analytic function of 2; in [/ and if there is a sequence of analytic 
functions aj{z) such that h{z;e) admits the series X^j>o as asymptotic expansion as e — » 

yz £ U and for all compact K C U there exists a constant C > such that |aj(z)| < C'''^^f Vz £ K. 
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Remark 2.3 We note that the function N{z) defined in (j2.19p is analytic in C\(— oo, 0] 
and it has poles on the negative axis. Using properties of the F-function, one can recover 
the limits 

lim N{z) = l, |argz|<7r, (2.20) 

|2|— >00 

lim N{z) = V2. (2.21) 

2— >0 

The poles of the reflection coefficient r(A; e) in a neighborhood of —1 correspond to the 

) . Since r(A) w — (A + 1) as A — > — 1, it is clear 

2 vre y 2 

that the poles of r(A;e) occur only for ImA < (cf. Remark 12. ip . although as e — > 
they are getting closer to —1. 

Remark 2.4 We observe that applying the condition (j2.7p we obtain for real A < 

|r(A;e)|2 + |t(A;6)|2 = l 
which shows that 

l-|r(A;e)|2 = e-2"(^)(l + 0(e)), as e ^ 0, (2.22) 

where r has been defined in (j2.15p . Therefore for —1 < A < 0, |r(A; e)p is equal to one 
modulo exponentially small terms. 

Let us now define 

K{X;e) = -r{X;e)ie^''^^\ as A G V+ := n+ U L'+(-l, 5). (2.23) 

where p has been defined in (j2.13p . 

It then follows that k(A) 1 uniformly in H^^. In D^{ — 1, 5), one obtains by (|2.18p 
that 

K{X;e) = ^7^(1 + 0(e)). 

For A G D_|_(— 1, 6), — iT(A) stays away from the negative real line, so that by (|2.2Up . 
k(A; e) is bounded as e ^ 0. 

Summarizing we can conclude from the above discussion that 

k(A; e) = l + 0{e), uniformly for A G n+ , (2.24) 

|K(A;e)| < M, as A G V+. (2.25) 

Example 2.5 Let us consider the potential uo{x) = —1/ cosh^x. In this particular 
case, the Schrddinger equation 



e'k. —f = A/, (2.26) 



1 

cosh'^ X 

has a solution of the following form, 



fix; A) = 2-'=(l - e)''/%Fiik -s,k + s + l,k + l;^il- 0), 



where 



^ = tanha;, k = ^^~2'^eV^~T' 
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and where sFi{a, b, c; z) is a hypergeometric function solving the equation z{l—z)w" + 
[c — {a + b+ l)z]'w' — abw = 0. As x ^ +00, we observe that ^ — > 1, so that 



f{x;X,e) ~ e" 



-Ax 



X +00. 



Transforming the hypergeometric function 

^ / , N r(c)r(c — a — b) ^ , , 

2Fi{a, b, c; z) = ) ' ^ -r:'^Fi{a, 6,a + 5 + l-c;l - 2;)+ 



V(c- a)V{c-b) 



r(c)r(a + fe - c) 

r(a)r(6) 

one obtains for x ^ —00 that 



(1 - zy-^-^Fiic -a,c-b,c+l-a-b;l- z), 



T{k + l)T{k) 



/lx,A,e; r(s + l)r(-s) ^ T{k-s)T{k + s + l) 

Therefore the reflection coefficient is equal to 
TC-V^ - s)TC-V^ + s + 



e e 



-Ax 



as X 



+ 00. 



r(A;e) 



(2.27) 



Using well-known asymptotic properties of the T -function it is a straightforward 
calculation to check that i2.12\} . {2.11), and \2.18\) hold for this example. 



3 Asymptotic analysis of the RH problem near the gra- 
dient catastrophe 

We will apply the Deift/Zhou steepest descent method on the RH problem for the 
KdV equation. The main lines of this method have been developed in [14^ I15j. but the 
analysis near the time and point of gradient catastrophe has not been established yet. 
The starting point of our analysis is the RH problem for M given in Section [2j 



RH problem for M 

(a) M(A; x, t, e) is analytic for A G 

(b) M satisfies the following jump conditions, 

/ 1 T-f \- f'lp2ia(A;a:,t)/e\ 

^^-(^)=^-WUa;.kU..)A T-Ux;e)\^ J' ^ < 

M+(A) = M_(A)c7i, ai = for A > 0, 



(c) M{X;x,t,€) 



for A 



00. 



-A 



-A, 



Our goal is to find asymptotics for M(A; x, t, e) in the small dispersion limit where 
e — > 0, while x, t tend at an appropriate rate to the point and time of gradient catas- 
trophe Xc, tc for the Hopf equation ut + 6uux = 0. 
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3.1 Construction of the ^-function 

A first crucial issue in the asymptotic analysis is the construction of an appropriate 
^-function. In [T3], a ^-function has been constructed when taking the small dispersion 
limit e ^ for fixed x and t. However we will, later on, take a double scaling limit 
where x and t tend to the critical values Xc and tc simultaneously with e — > 0. For 
this purpose we need to modify the ^-function in such a way that we will be able to 
construct a local parametrix for A near the critical point Uc- 

Let Q be defined as follows, with Uc = u{xc,tc), 

g{\;x,t) = / , di], (3.1 

where p and a are defined as before by 

a-_(A) 



p{X) = x-{X)\^+ iVuoix) - A - V^]dx, 

a{X; x, t) = 4t{-Xf/^ + x{-X)^^\ 



and \/uc — X is analytic off [uc+oo), positive for A < Uc- Now Q is analytic in C \ 
[uc, +oo) and ^(A) = C(A~^/^) as A ^ oo. One checks directly that 



gi{x,t):= \\m{-XY''g{X-x,t) = - ' ' -dr]. (3.2) 

Since p{rj) is independent of x, we obtain 



dxQi{x,t) = — [ J —dr] 

_ Uc 

2 ■ 



(3.3) 



In the following proposition, we collect two boundary value relations for the Q- 
function. Together with the asymptotic behavior of Q, we could have considered those 
as the defining properties for the Q as well, instead of defining it explicitly by means of 
(j3.ip and deriving the properties afterwards. 



Proposition 3.1 G is analytic inC \ [uc+oo), and we have 

g+(A)+g_(A) = 0, /or A G (0,+oo), 

g+(A) + g_(A) - 2p{X) + 2a{X) = 0, for X G {uc, 0), 

g+(A)-g_(A) = 0, for Xe i-oo, Uc). 

Proof. The first statement is immediate since \/uc — Xj^ = —\/uc — X_ and since 
A > does not belong to the interval of integration. For the second statement, note 
that as A S (iic,0), ^+ + Q- is equal to a contour integral around A. A simple residue 
argument gives us the required result. The last statement is obvious because \Juc — X 
is branched for X> Uc and X <Uc does not belong to the interval of integration. □ 

3.2 First transformation of the RH problem 

In this section we perform a first transformation of the RH problem for M. The goal 
of this transformation is to simplify the jump matrices for the RH problem, using the 
properties of the ^-function in an appropriate way. 
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We define 

T{X;x,t,e) = M(A;2:,t,e)e-^^(^'^'*)'"^ 



(3.4) 



From the RH conditions for M and the definition of T, we can now derive the RH 
conditions for T. 

RH problem for T 

(a) T is analytic in C \ M, 

(b) T+(A) = T^{X)vt{\), as a G M, 

VTiX) = et^?-W'^3^;A/(A)e-^^+W'^^ 

(c) r(A) ~ as A ^ oo. 

Using Proposition 13.11 one checks that the jump matrix vt transforms to 



vt{X) 



0-1, 



r{X)e' 



-f(A)e-f 1 - |r(A)|^ 



as A > 0, 

as A G (— oo, Uc), 



e E 



i(g+(A)-g_{A)) 



r(A)e' 



V{A) 



-r{X)e-fpW (1 - |r(A)|2)et(^+W-^-W) 



as A G {uc, 0), 



(3.5) 



It is important that we can express the KdV solution u{x, t, e) in terms of the new RH 
problem for T. Using (j2.1ip . (j3.3p . and (j3.4p . the solution of the KdV equation can 
now be obtained from 



u{x,t,e) = 2dxQi{x,t) -2iedxTl^{x,t,e) 
= Uc- 2ied^Tli{x,t,e), 

where T^i is given by 



(3.6) 



Tii(A;x,t, e) = 1 + 



+ 0(A-i) 



as A — > oo. 



Let us define an auxiliary function </> in V = V+ U V_ , where V+ is the region defined 
in (12:23]) and (l2l6D . and V_ = V^, by 



(j){X; X, t) = G{X; x, t) - p{X) + q(A; x, t), 



(/)(A; X, t) = (/>*(A; x, t) = (/)(A; x, t), 



for A G V+, 
for A G V-. 



(3.7) 
(3.8) 



Note that because of reality, is analytic across (— l,Uc), but not on (—1 — (5,-1), 
because this is a part of the branch cut for p, as we discussed in Section [2j On (uc, 0), 
(p is not analytic because of the branch cut for Q. We obtain from Proposition 13. II that 
the following equation holds, 



2(/)_|_(A; X, i) = Qj^{X;x,t) — Q-{X;x,t), for A G (uc,0). 



(3.9) 
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Lemma 3.2 The function (p defined in J^. 7^-[J78\) takes the form 



cP{X; X, t) = - A(x -Xc- Gucit - tc)) + / {fUO + Qt)^/C - Xd^. (3.10) 



Proof. Using the definition (j3.1|) of Q and (j2.13p . it follows that 



vr 



(r/ - \)yjr] - Uc 



drj. 



For A ^ (tic, 0) one uses the residue theorem to integrate the terms in x and t while for 
the double integral we exchange the order of integration obtaining 

g{\) = yjuc - \{x - Quct) + At{uc - A)^/^ - a(A; x, t) 



+ 



■y/ttc - A 
27r 



drj 



:d^. 



Applying the residue theorem to the last integral one obtains 



^(A) = v^Uc - X{x - Quct) + 4t{uc - A)^/^ - a(A; x, t) 



+ 



+ + 



Integrating by parts the last term of the l.h.s. of the above relation and using the 
identity Xc = Quctc + f-{uc), one arrives at the statement. In the case A G {uc,0) one 
has to evaluate the principal value of the integral defining Q in the same way as done 
above. □ 

The introduction of enables us to rewrite the jump matrix vt in a more convenient 
form. 

Indeed we have that 

-71, 



1 r(A)e^(^(^)+"(^)) 

2i 



1 



iK+(A)eT<^+W 



iK*_(A)e-T^-W l-|r(A)|^ 



e e 



0+{A) 



iK{X) 



\ iR{X) (1- |r(A)|2)eT<^+W^ 
As before we have written k for 

k(A; e) = -ir{X; e)e^''^^\ as A G V4 



as A > 0, 
as A < —1 — 5, 

as A € (—1 — 5, Uc 

as A S {uc, 0). 



(3.11) 



(3.12) 



with boundary values on M denoted by k+. We have written k'L{X) = R-^-{X). In 
particular near Uc, where our main focus is, the introduction of (f) will turn out to be 
convenient. 



In the following proposition, we discuss the behavior of (p in different regions of V. 
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Proposition 3.3 For sufficiently small (5 > 0, there exists (5i > and there exists a 
neighborhood of {—1 — 6,Uc — 6) such that for \x — Xc\ < 5i, \t — tc\ <Si, the following 
holds, 



(3.13) 



lm(j){X;x,t) > 0, as X £W r\{z £ C -Amz > 0}, 

Im (j){X; x,t) <0, as X £W n{z £ C -.Imz < -0}, 

Im 0+(A; X, t) < 0, as X £ [uc + 6,0], 

-r(A) + i(j)+{X;x,t) < 0, as A G [uc + 6,0]. 

Proof. Taking the derivative of (j3.10p one obtains 

(/.'(A; x, t) = - Jl J x -xc- Qucit - tc)) - \ T d^. (3.14) 

Since /i(C) + Qtc < for ^ G (-1,0), it follows that 

(l)'{X;Xc,tc) > 0, for X £ [-l,Uc- 6]. 
By continuity there exists 6i > such that for |x — xd < 6i, \t — tc\ < 6i, 

(l)'{X;x,t) > 0, as X£[-l,Uc-6]. (3.15) 

From the above inequality, by the Cauchy-Riemann relations, it follows that 

Im0(A;x,t) > 0, for A G W n {A G C : Re A > -1, Im A > 0}, (3.16) 

lm(j){X;x,t) <0, for A G Wn {A G C : ReA > -1,1mA < 0}, (3.17) 

at least if the neighborhood W is chosen sufficiently small. 

In the case A < — 1 close to —1, we need to be a little bit more careful because of 
the branch cut for (p. The function /_ , the inverse of the decreasing part of the initial 
data, now assumes the form 

/_(A) = XM - -^^={1 + 0{X + 1)), as A ^ -1. 

Since the function \/A + 1 is analytic in C\(— cxd, —1] and positive for A > —1, it follows 
from (f3T0]) that 

Im0+(A;x,t) > 0, for A G [-1 - 5, -1), (3.18) 

Im0_(A;x,t) < 0, for A G [-1 - (5, -1), (3.19) 

Since Re i?^'-|-(A; x,t) > for A G [—1 — 5,-1], by the Cauchy-Riemann relations it again 
follows that the first and second inequalities in (j3.13p hold true. 

Regarding the third inequality in (j3.13p we have for A G [uc + 6, 0] that 
i(/)+(A; X, t) = ^/X- uc{x -Xc- ducit - t,)) - / (/i(0 + 6t) y^A^de. 

J Uc 

This means that i(p+{X; Xc, tc) > for A G [uc + S,0]. When \x — Xc\ < 5i, \t — tc\ < Si 
with 6i > sufficiently small, it follows that i(p+{X; x,t) > for A G [uc + 6,0] . 
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For proving the last inequality of Proposition 13.31 with r defined as in (j2.15p . it is 
straightforward to verify that for A G [«c + S, 0] , 

— t(A) + icl)^{\; X, t) =\/ \ — Uc[x — Xc — 6uc{t — tc)] — 4t(A — Uc)^ 

where /+ is the inverse function of the increasing part of the initial data uo{x). Since 

—t{X; Xc,tc) + {X',Xc, tc) < 0, as A G [uc + 6,0], 

there exists 6i > sufficiently small, such that —T{X;x,t) + i(j)^{X; x,t) < when 
\x — Xc\ < Si and \t — tc\ < Si. □ 

As e ^ 0, the above proposition shows us exponential decay of the diagonal entries of 
vt on [uc + f^, 0], and oscillating behavior of the off-diagonal entries oi vt on {—1 — 6, Uc) ■ 
Performing a next transformation will enable us to deform the oscillatory entries to 
exponentially decaying entries as well. 



3.3 Opening of the lens 

We are able to factorize the jump matrix in the following way for —1 — S < X < Uc 



vt{X) 



m!L(A)e 




iAc+(A)eT<^+W 



1 



(3.20) 



Furthermore the first factor can be extended analytically to the lower half plane and 
the second factor can be extended to the upper half plane. This observation enables us 
to move the jump contour, which coincided with the real line so far, into the complex 
plane. The spirit of the Deift/Zhou steepest descent method is that, deforming the 
contours, one can deform oscillatory jump matrices (on the real line) to exponentially 
decaying jump matrices (in the complex plane). 

We will use the factorization of the jump matrix to open lenses along the interval 
(— 1 — S,Uc) for some sufficiently small but fixed S > 0. It is not necessary to open 
lenses elsewhere on the real line, because there the jump matrices will turn out to be 
exponentially small already without deforming the contour. However it is necessary 
to open the lens starting from some point slightly to the left of —1. Opening the lens 
exactly at —1 would not lead to exponentially small jump matrices near —1. Let us 
consider a lens-shaped region as shown in Figure [TJ We retain the freedom to specify 
the precise choice of the lens later on, bur for now we assume that Si is, except near 
Uc, contained in W n {z G C : Im z > 0}, and that T,2 is, except near Uc, contained in 
W f] {z £ C : Imz < 0}, where W is a region for which Proposition 13.31 is valid. 

Define S as follows, 



S{X) 



T(A) 



-in{X)e 
1 



in region I, 



( 0(A) ij' in region II, 



(3.21) 



nx). 



elsewhere. 



with K*{X) = k{X). 

Now the RH problem for S takes the following form. 
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Figure 1: The jump contour S5 and the jumps for S 



RH problem for S 

(a) S is analytic in C \ Ss, 

(b) S+{X) = S^{X)vs for A G Es, with 



1 0' 



(c) Six) 



1 



\ IK 

vt{X), 

1 



</'+(A) 



iK,{X) 



(A) (1 - |r(A)|2)ef i ' 



on Si, 
on T,2, 
as A S (uc, 0), 



(3.22) 



as A G (-00, -I - 5)U (0, +00). 



-A 



as A ^ 00. 



Since S{X) = T(A) for large A, formula 
the KdV solution by 

u{x, t, e) = Uc- 2iedxSli{x, t, e), 



remains valid for S, so that we can retrieve 



(3.23) 



where 



Sii(A;x,t,e) = l + ^Sli^^ + 0(A-i), as A ^ 00. 



(3.24) 



In the following proposition, we show that the jump matrix vs converges to a 
constant matrix uniformly fast as e ^ 0, except for in a neighborhood of Uc- We should 
note that this is only true because of the semiclassical limit of the reflection coefficient 
described in Section [2l so indirectly this relies on Assumptions [TTTl 

Proposition 3.4 We can choose a suitable contour S5 such that for any neighborhood 

U of Uc, 



i;s(A) =t;(°^)(A)(I + 0(e)), as e ^ 0, 
uniformly for X G S5 \ U, with v^°°^ defined by 



?;(°°)(A) 



fJi, forX>0, 
iai, for X £ {uc,0), 
/, elsewhere. 



(3.25) 



(3.26) 
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Proof. (i) For A > 0, ^^(A) = v^°°\X) so that the result holds trivially here. 

(ii) For A G {uc, 0) \U, it follows from (j2.24p that that the off-diagonal entries of vs{X) 
tend to i with uniformly small error of order e. The 11-entry and the 22-entry of 
vs{X) are exponentially small because of (j2.15p and Proposition 13.31 

(iii) For A G (Si U S2) \ U, under the assumption of a well chosen contour in view of 
Proposition 13.31 it follows from ()2.25p and Proposition 13.31 that 

^;5(A) =t;(°°)(A)(/ + 0(e-i')), ase^O, 

for some constant C3 > 0. 

(iv) For A < —1 — 5, vs{X) I because of the uniform convergence of the reflection 
coefficient, see (|2.17p . together with the reality of Q and a. 

□ 

Remark 3.5 Note that the error term of the jump matrices is exponentially small 
everywhere except for on the interval {uc,0), where the WKB approximation for the 
reflection coefficient causes the C'(e)-error. 

3.4 Outside parametrix 

Ignoring the exponentially small jumps and a small neighborhood U of u^. where the 
uniform exponential decay of the jump matrices does not remain valid, our RH problem 
reduces to the following RH problem for 

RH problem for P(°°) 

(a) P(°°) : C \ [uc, +00) ^ C2x2 is analytic, 

(b) p(°°) satisfies the following jump conditions on (uc, +00), 

p{oc) ^ p^oo)^^^ ^ ^ ^^^^ ^3_27) 

pjoo) ^ ipi°°)^^^ as A G {uc,0), (3.28) 

(c) p(°°) has the following behavior as A — > 00, 

We can construct the solution of this RH problem explicitly as follows, 

p(-)(A) = (-A)i/^(u,-A)-'^3/4(^l (330) 

from which it follows directly that the asymptotic condition (|3.29p can be strengthened 
to 

pM(A) = (/ + ^cT3 + 0(A-2))(.^_^^i/2 _i(_\)i/2), as A -00. (3.31) 

In the outside region away from Uc, the leading order asymptotics of S will be 
determined by p(°°). To obtain uniform asymptotics and asymptotics beyond the 
leading term, we still need to construct a local parametrix near Uc which matches with 
the outside parametrix at dU. This is the goal of the next section. Note also that the 
uniform decay of the jump matrices remains valid near —1 and 0, so that there is no 
necessity to construct local parametrices near those points. 
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3.5 Local parametrix near Uc 

The uniform convergence of the jump matrices as e ^ breaks down near u^. Here we 
need to construct a local parametrix by mapping a suitable model RH problem onto 
a neighborhood lA of Uc- We will use a model RH problem which is associated with 
a fourth order analogue of the Painleve I equation, the equation (jl.Sp . The aim 
of this section is to construct a local parametrix P in U which has approximately the 
same jumps as S has in U, and which 'matches' with at dU. With 'matching' 

we mean that P{X)P^°°\X)-^ tends to the identity matrix in a suitable double scaling 
limit. The appropriate double scaling limit will turn out to be the one where we let 
e ^ and in the same time we let x ^ Xc and t — > tc in such a way that t — tc = C'(e^/'^) 
and X — Xc — 6uc{t — t^) = 0{e^^'^). More precisely, we want P to satisfy a RH problem 
of the following form. 



RH problem for P 

(a) P : Z7 \ ^ C2><2 is analytic, 

(b) P satisfies the following jump condition onU r\ S5, 

P+(A) = P_(A)^;p(A), 
with vp given by 



vp{X) 







le 



1 



-f</.(A;x,t) ]_ 



^g-f 0+(A;x,i) • 




as A G {uc, 0), 



(3.32) 



(3.33) 



(c) if we perform the double scaling limit where we let e — > and at the same time 
we let X — > Xc and i — > ic in such a way that 

6(t - tc) 



X — Xr — QUrit — tr) 

lim — = X, 



lim ■ 



(8A:e6)i/7 ' (4A;3e4)i/7 

the following matching condition holds, 

P(A)P(°°)(A)-1 ^ /, forAeaZ^. 

We will show later on that vp approximates vs in the following sense, 

vs{z)vp^{z) = I + 0{e), uniformly for z e n as e ^ 0. (3.36) 



(3.34) 
(3.35) 



3.5.1 Model RH problem associated with Pf 

Recall that the Pf equation is the following fourth order differential equation for U 
U{X,T), 



X = TU 



1. 



1 



1 



(3.37) 



It was conjectured in [_17] and proven in p] that this equation has a real pole- free solution 
U {X, T). In agreement with the conjecture of Dubrovin, it follows from previous results 
in [32\ [33] that the pole-free solution is unique. The RH problem characterizing this 
solution is the following, see [271 E] • 
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III 



Figure 2: The jump contour for P 

RH problem for ^: 

(a) ^ = ^'(C; X, t) is analytic for (" G C \ F, with F as shown in Figure O 

(b) ^ satisfies the following jump relations on F, 

M/+(C) = ^-(C)(_°i J), forCGTs, (3.38) 

*+(C) = ^-(C)(j I), forCGTi, (3.39) 

^+(C) = ^-(C)(J J), forCGF2UF4. (3.40) 

(c) ^ has the following behavior at infinity, uniformly for {X, T) in compact subsets 
of \ V, where V denotes the set of poles of U , 

*(C) = r4-A'(/+«.3r"'^ + i(_'^; 5)<- 

+Ra3C'l^- + 0(e)] e'"!*'-'^''"!", (3.41) 



where U = U{X,T) is the real pole- free solution of the Pf equation (j3.37p . 
dxQ{X,T) = U{X,T), R is some unimportant function of X and T, and 
and 9 are given by 



N 



i= (_\ I) e-x-, eiC;X,T) = r^.C'^' -^C'^' + X^^^ (3.42) 



Remark 3.6 It is important for us that, for each {X,T) G R^, there is a neighborhood 
^ C of {X, T) such that ^'(.; x, t) is well-defined for all {x^t) ^ A and such that the 
asymptotic condition (j3.41|) holds uniformly for (rr, t) S A. This is only guaranteed by 
the fact that U {X, T) has no real poles [8] . 

Remark 3.7 There is some freedom in choosing the angles that F2 and F4 make with 
the negative real axis, as long as those angles are in absolute value smaller than ^ 
(which is the angle for the Stokes lines closest to the negative real line). For simplicity, 
we fix F2 = {2; € C : argz = ^} and F4 = {z G C : argz = ^|^} to be the so-called 
anti-Stokes lines. 
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Figure 3: The jump contour for T for ^ 



The RH solution ^ provides in each sector fundamental bases of solutions to the 
Lax equations associated to the Pf equation. Those equations are 

= A-if, -ifx = B-^, = C*, (3.43) 

where the matrices A, B, and C are polynomials in given by 

f-^UxC - {l2UUx + Uxxx) Se + sue + + 2Uxx - 120T) 

~ 240 A21 AUxC + {l2UUx + Uxxx) 

A21 = 8C^ - 8UC^ - (4C/2 + 2Uxx + 120T)C + (IGU^ - 2U^ + AUUxx + 2AQX), 

/ 1\ l( Ux -2C-2U\ 

\C-2U Q)' Q\-2C^ + 2UC + c{X,T) -Ux )' 

Compatibility of the first and the second equation in (|3.43p shows (after a not so short 
derivation) that U solves the P]' equation, while compatibility of the first and the third 
equation in (|3.43|) shows that U{X,T) solves the KdV equation normalized to the form 
Ut + UUx + = 0, see also [IT]. One can make the rescalings as in formula 

(II. 7p to obtain the KdV equation in the form (II. 3p . 

3.5.2 Modified model RH problem 

In order to obtain a RH problem that we can use to construct the local parametrix 
near Uc, we transform the RH problem for ^ by defining <I> in the following way, 

feX-3vi/(^)ee(C;X,T)<x3 asImC>0, 
mX,T) = l V V (3.44) 

[e-x'^3^((-)g9(C;^,T)a3g2iia3^ as ImC < 0. 

Clearly <I> is analytic in C \ F, where F is the jump contour for ^ shown in Figure [3l 
Using the fact that 9+{C;X,T) = -$-{(■, X,T) for ( e (-00, 0) and dOS]) . one can 
check using the definition (13.441) that ^ has no jump on (— oo,0) and is thus analytic 
on (—00, 0). Using (j3.39p . (j3.40p . and (|3.44p . one checks that the jumps along the other 
rays in the contour do not vanish but are modified by the transformation. By (j3.4ip 
and (j3.44p . also the asymptotic behavior of ^ is different from the asymptotic behavior 
of ^. $ satisfies the following RH problem. 
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RH problem for $ 

(a) $ is analytic for C G C \ f , with f = Ti U r2 n 

(b) $ satisfies the following jump relations on T, 

a>+(C) = 1>-(C)(^ . qJ, forCeri, (3.45) 

$+(C) = '3>-(C)(^Q J, forcers. (3.46) 

$+(C) = $-(C)(.g2e(J;X,T) forCer4. (3.47) 

(c) $ has the following behavior at infinity, 

$(C) = e-x'^s^-i^s^ 

e"*^^, as C ^ oo with Im^ > 0, 



X < 



-1 

1 



(3.48) 



e 4 as C — > oo with Im^ < 0, 



which can be rewritten in the following way, both for C — > oo in the upper and 
the lower half plane, 

.,c,.-L,-o-j"(_\ ;)(...«..(-c)--i('^ 

+tRa3{-Cr^^^ + 0{C^)) . (3.49) 



This behavior holds uniformly for {X,T) in compact subsets of \ V, where V 
is the set of poles of U. 

3.5.3 Construction of the parametrix 

We are now ready to specify the form of the parametrix P. Let P be of the following 
form, 

P(A) = E{X; e)cl>(e-2/V(A); e-%i(A; x), e'^" g2{\; t)), (3.50) 

where E, /, gi, and g2 are analytic in U and x = x — 6uct. Furthermore we require that 
/(M nU) CM with /(uc) = and f'{uc) > 0, in other words / is a conformal mapping 
from ^ to a neighborhood of the origin. We will determine the precise form of E, /, 
gi, and g2 later. 

First we specify our choice of the contour S5 n^/ by requiring that /(S5 ClU) C F, 
which we can do since / is a conformal mapping. By this construction P and S have 
their jumps on the same contour in U. We will now use the freedom we still have in 
defining f,gi, and g2 in order to create jumps for P that are the same as the ones 
specified in (j3.32p - (|3.33p . Afterwards we will define E in such a way that the matching 
condition (j3.35p is satisfied as well. 
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3.5.4 Definition of /, gi, and g2 

Our goal is to define /, gi, and (72 in such a way that 

^(e-2/7/(A); e-6/V(A; x), e'^/^^alA; t)) = i0+(A; x, t), for A G (n^, 0) n U, 

0(e-2A/(A);e^6/7^i(A;x),e-4/752(A;t)) = i0(A;5,t), for Xe^iDU, (3.51) 

0(e-2/7/(A); e-6/V(A; x),e-^/-' g2iX; t)) = -\^{\- x, t), for A G S2 n 

with (/) and defined by (j3.7p and (j3.42p . For simphcity, we changed variables for 
and consider it now as a function of x instead of x, with 

X = X — 6uct. 

We should now compare the jump matrices (j3.33p for P to the 'model' jump matrices 
(|3.45p - (|3.47p for <I>. The above conditions (|3.5ip imply that the parametrix P, defined 
as in (j3.50p . satisfies the requested jump condition (j3.32p . 

With 6(0 having its branch cut along the negative real line, one verifies that the 
equations in (|3.5ip are satisfied if 

^(-/(A);-5i(A;S),52(A;t)) = -0(A;x,t), for A G Z^, (3.52) 

In agreement with this condition, we first define / in such a way that 

Y^(-/(A))'/' = -0(A;xe,tc). (3.53) 
Using (11. 9p and (13.101) . integrating by parts twice leads to 



(/)(A; X, t) = \Juc - X{x - Xc) + 4(mc - A) 2 (t - tc) 

4 

+ 15 



f^^ f'-{m->^)^di. (3.54) 



Keeping track of the branch cuts for the fractional powers, this implies that / is analytic 
in lA with 

/(n,)=0, /'(^xj = (8A;)?, k = -f!!{uc)>Q. (3.55) 
Now we can define 32 by requiring that 

^92(A;t)(-/(A))3/2 = 0(A;Xe,t) - 0(A;ie,tc) = ^{uc - A)t(t - Q. (3.56) 
Since / is a conformal mapping, this defines g2 analytically in ZY, with 

52K;t) = |^. (3.57) 

Finally we define gi by the equation 

gi{\; X, t)(-/(A))i/2 = ^(A; i, t) - (/)(A; ie, t) = ^ Uc - A(x - x^). (3.58) 

Again using the fact that / is a conformal mapping, g\ is analytic in lA^ with 

I ~ ,^ X - Xc - 6uc(t - tc) ,„ 

5ri(uc;x,t) = J . (3.59) 
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Summing up (|3.53|) . (j3.56p . and (j3.58p . we find using (|3.42|) that indeed condition (j3.52p 
is satisfied. Tliis means tliat the jump conditions (|3.32p - ()3.33p for P are valid. 

The function <I>((^;X, T) is not defined for values of X and T where U{X,T) has a 
pole. In order to ensure that the parametrix is well-defined, we need to know that U 
has no pole at {e^^/'^gi{X;x,t),€~'^^'^g2{)^;t)). In the double scaling limit where e ^ 
and at the same time x — > Xc and t ^ tc in such a way that 

X - Xc-6uc{t-tc) ^ Q{t-tc) „ 

hm , . ^, , ,„ = A , lim ■ — 



(8A;e6)i/7 ' (4A;3e4)i/7 

it follows from (|336D - (l339]) that 

e-6/V(A;x,t) ^X, e-4/7^2(A;t) ^T, as A ^ u,. 

Prom the fact that U{X, T) is meromorphic both in X and T and that it has no poles 
for real values of X, T, it follows that there is a pole-free neighborhood of {X, T) in 
(cf. Remark [3.6p . in which (e~®/^gi(A; x, t), e'^^"^ g2{^\ t)) is contained for sufficiently 
small e, provided that A lies in a sufficiently small neighborhood U of Uc- 

In order to satisfy also the matching condition (j3.35p . we define the analytic pre- 
f act or E by 

E{\; e) = ^P'^W (I (-6-2/V(A))^. (3.60) 

Using the definition (j3.30p of P^°°\ one checks directly that E is analytic in U, as it 
should be in order to have a suitable parametrix. The matching condition (j3.35p for 
A € dU in the double scaling limit can now be specified as follows, 



P(A)P(°^)(A)-1 = P°°(A) (/ + iQcT3(-/(A))-V2ei/7 

(3.61) 

Here we used the abbreviations 



U = C/(e-6/V(A; x, t), e-^/'g2{X; t)), (3.62) 

Q = Q(e-6/V(A; x, t), e-"/V(A; t)), (3.63) 

R = i?(e-6/V(A; X, t), e-^/^52(A; t)). (3.64) 

This ends the construction of the local parametrix. 

3.6 Final transformation 

We define 

|5(A)P<».(A)-, ^X,C\U. 

^ ^ \S{X)PiX)-\ asXGU. ^ ^ 

Note first that outside parametrix has been constructed in such a way that R has no 
jump on (0, +00). For z £U H S5, we have that 

RZ\z)R+{z) = P^{z)vsiz)vp\z)Pz\z). (3.66) 
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Figure 4: The contour S/j after the third and final transformation. 



On one hand it fohows from the construction of the parametrix that P_ (z) is uniformly 
bounded for z £U D S5. On the other hand 



vsvp 



1 i{K — l)e" 
,0 1 



1 

l)e- 

K 



i{K — l)e 
+ (1 - 



on Si n U, 
on S2 n Z//, 
on (ucO) nZ^. 



(3.67) 



Except for the 21-entry on (uc, 0), the exponentials in the above matrices are uniformly 
bounded on the jump contours inside U because of (|3.5ip and the fact that the jumps 
for <I> are uniformly bounded on the jump contour T. Here the WKB approximation 
(|2.24|) for K ensures that vsVp^ = I + 0{e). If in addition we use (j2.22p and Proposition 
13.31 the same follows for the 21-entry on (uc,0) CiU. By Proposition 13. 4^ one checks 
that the RH problem for R has the following form. 

RH problem for R 

(a) R is analytic in C \ E/j, with E/j = (S^ U dU) \ (0, +00) as shown in Figure HI 

(b) R+{X) = R-{X)vr{X), where the jump matrix vr has the following double scaling 
asymptotics, 



vr{X) 



'P(A)P(°°)(A)-1, asXGdU, 
/ + 0(e), as A G \ (9^/. 



(3.68) 



(c) R{X) I as X ^ 00. 



This is the final RH problem in our analysis. In the next section, we will obtain 
uniform asymptotics for R which will enable us to prove Theorem 11.21 



4 Proof of Theorem 11.21 



Using (|3.61|) . we can expand the jump matrix vr in fractional powers of e in the double 
scaling limit, 



vr{X) = 1 + 6i/^A(i)(A) + e2AA(2)(A) + 0{e'/'), 



(4.1) 
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with 

A«(A) = iQ ■ (-/(A))-i/2p{oo)(^y^p(oo)(^)-i^ ^4 2) 

A(')(A) = -^(-/(A))-ip(-)(A)(^' ^^P^^\X)-\ (4.3) 
for A G dlA^ and 

A(i) (A) = A(2) (A) = 0, for A G S/j \ (4.4) 

since the jump matrices are equal to / up to an error of 0{e) on the other parts of the 
contour. Note that the functions A^-*^) and A^^^ are meromorphic functions in U with 
poles at Uc- 

It is a well-known result that the uniform asymptotic expansion (14. ip for the jump 
matrix yields an asymptotic expansion of the same form for the RH solution i?, uni- 
formly for A G C \ Sr, 

R{X) = 1 + €^'^R^^\\) + e^/^R^'^\X) + 0(e3/7). (4.5) 

This can proven exactly as in [111 [13] by estimating Cauchy-type operators associated 
to the RH problem. 

Remark 4.1 We should note that, in [HI [121 US]) the local parametrices erase the 
jumps inside disks near the special points completely. Because of the 0(e)-error of the 
reflection coefficient to its WKB approximation, we have not been able to remove the 
jumps in lA completely. This leads to a final RH problem which is similar to the one 
in [2]. However the remaining jumps in the interior olli are of smaller order than the 
relevant jumps on dU, and therefore this does not cause any further problems. 

The compatibility of the expansion (j4.1l) for the jump matrix with the expansion 
()4.5p for the RH solution leads, using the jump relation R+{X) = R-.{X)vfi{X), to the 
following conditions for A G dU, 

(A) = rW (A) + a(i) (A) , (4.6) 
i?f (A) = R^^\X) + r'l\x)A^^\X) + A(2)(A). (4.7) 

Note in addition that i?(A) ^ / as A ^ oo, and thus also R^^\X) — * / for j = 1, 2. The 
jump and asymptotic conditions constitute additive RH problems for i?^^^ and R^'^\ 
By inspection we see that their unique solutions are given by 

^a)(,) f^Res(A(|);..), asAGC\^ 
\^Res(A«;u,) - AW(A), as A G Z^, 

i?(2)(A) = 

r_J_Res(ii(i)A«+A(2);n,), asAGC\Z^, 
|_J_Res(i?(i)A«+A(2);n,)-i?«A«(A)-A{2)(A), as A G Z^. 

After a straightforward calculation we find using (j4.2p . (j4.3p . and (j3.30p that, for A G 
C\Z^, 

= -Q/>.r'/^^(» (4.10) 
««(A) = (: »). (4.11) 
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where the *'s denote unimportant entries, and where we have now written 

We observe aheady that the main sub-leading terms in the asymptotic expansion for R 
are determined completely by the matching of the local parametrix P with the outside 
parametrix 

Besides the expansion in negative powers of e, R admits also an expansion in neg- 
ative powers of A as A — > oo jTT], 

R(X) = I + ^ + 0{X-^), asA^oo. (4.12) 
A 

Compatibility of the small e-expansion (14. 5p with the large A-expansion ()4.12p learns 
us the following in the double scaling limit, 

i?i,i2(x,t,e) = -6i/7Q/'K)-i/2 + 0(63/7). (413) 

Now it turns out that the KdV solution u{x,t,e) is contained in Ri. Indeed by (13.651) . 
we have that 5(A) = i?(A)P(°°)(A) for large A, which means by (f3:§T]) and KT2\\ that 

^iiW = l + ^(T^ + ^(^~')' asA^oo. (4.14) 

In view of (j3.23p we obtain using dxQ = U that 
u{x,t,e) = Uc- 2iedxRi,i2 (5, t, e) 

/o.2\ 1/7 

= (^J [/(e-6/V(«c;x,t),e-'/'52(ne;t)) + 0(e^/^). 

By p.57p and (j3.59p , this completes the proof of Theorem 11.21 
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